Electron fractionalization and unconventional order parameters of the t-J model 
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In the t-J model, the electron fractionalization is unique due to the non-perturbative phase string 
effect. We formulated a lattice field theory taking this effect into full account. Basing on this field 
theory, we introduced a pair of Wilson loops which constitute a complete set of order parameters 
determining the phase diagram in the underdoped regime. We also established a general composition 
rule for electric transport expressing the electric conductivity in terms of the spinon and the holon 
conductivities. The general theory is applied to studies of the quantum phase diagram. We found 
that the antiferromagnetic and the superconducting phases are dual: in the former, holons are 
confined while spinons are deconfined, and vice versa in the latter. These two phases are separated 
by a novel phase, the so-called Bose-insulating phase, where both holons and spinons are deconfined 
and the system is electrically insulating. 

PACS numbers: 74.40.Kb, 74.42.-h 

I. INTRODUCTION 

One of the canonical notions in the theory of thermal phase transition is the Landau-Ginzburg- Wilson (LGW) 
symmetry-breaking paradigm. For a second order phase transition, one utilizes a few "order parameters" to charac- 
terize different symmetries of the phases and thus a critical theory is in terms of fluctuations of these order parameters, 
which drive the system towards a critical point. In the past few years, whether the LGW paradigm may be extended 
to quantum phase transitions has become one of the most fascinating topics in studies of strongly correlated electronic 
materials. In fact, there has been substantial investigations (e.g., Refs. [TH3] ) showing that the conventional LGW 
paradigm fails to describe diverse quantum phase transitions. In particular, a microscopic mechanism that violates 
the LGW paradigm is the so-called "fractionalization", discovered, say, in the study of two-dimensional quantum 
magnets [1] . It suggests a very general scenario of quantum phase transitions that has been currently under intensive 
experimental and theoretical investigations. Basically, a quantum particle may be effectively fractionalized into a few 
degrees of freedom interacting via emergent gauge field, and some fractionalized degrees of freedom appear at the 
critical point but not in the stable phases at each side of the critical point. 

Thus, as correlated electronic systems are concerned, a fundamental issue is: how does the electron fractionalization, 
where electrons are effectively fractionalized into the spin and the charge degrees of freedom at short scales, lead to 
various phases and the phase transitions between them? In one dimension, it has been well established that the 
spin-charge separation results in the Tomonaga-Luttinger liquid. In contrast, less has been known about the interplay 
between the electron fractionalization and quantum phase transitions, notably in the two-dimensional doped Mott 
insulator [5] which is widely accepted to be the prototype of high temperature superconductors [6-8 J . Notwithstanding, 
experimental and theoretical analysis indicating the electron fractionalization to be at the root of low-dimensional 
strongly correlated electronic systems has been well documented. In fact, there exist two classes of experiments 
on cuprate materials, namely the large Nernst effect in thermoelectric transport measurements [2j and the periodic 
modulations of electronic density of states in scanning tunneling microscopy studies [10] . It has been pointed out [3] 
that a LGW-type theory trying to reconcile these two phenomena cannot be, at least, extended to zero temperature. 
Instead, substantial theoretical analysis shows that the electron fractionalization necessarily leads to new quantum 
phases 6 . This opinion has been reiterated very recently by Zaanen and Overbosch [TT], who pointed out that the 
drastic change in the nature of quantum statistics upon doping - a direct manifestation of the electron fractionalization 
- may be the key to understanding physics of the pseudogap phase. Central to these studies are a number of 
fundamental issues. The most prominent ones include: What kind of quantum statistics do the emergent degrees of 
freedom obey? How do they interact with each other? How does the electron fractionalization turn an antiferromagnet 
into a superconductor upon doping? What does the phase diagram look like? 

While to investigate these issues for general doped Mott insulators proves to be highly challenging, in this work, we 
switch to a limiting case where the on-site Coulomb repulsion is strong enough, namely the t-J model (on a bipartite 
lattice), and will present an analytical study of the above problems. Our starting point is an exact non-perturbative 
result for this model, the so-called phase string effect discovered some time ago |12) . Specifically, the motion of a hole 
may leave a trace on the lattice plane upon colliding with spins. The phase string effect means that such two nearby 
propagating paths of the hole, though possessing (almost) identical amplitudes, may substantially differ in the phases 
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by 7r (namely in the overall sign) . Such a sign structure accounts for the parity of the number of hole-spin collisions 
and is thereby intrinsic to individual paths. Most importantly, with this effect being taken into full account, it turns 
out (to be detailed in this work) that an electron is uniquely, rather than arbitrarily, fractionalized into two bosonic 
constituents namely the holon and the spinon, satisfying a mutual statistical interaction. Basing on this exact result, 
in this work, we will study-at a full microscopic level-how a phase diagram flows out of electron fractionalization. 

The main results of this work were briefly reported in our earlier publication |13j . The present paper is devoted 
to reporting the technical details and the remaining is organized as follows. For the self-contained purposes we will 
briefly review the phase string effect in Sec. [IT] and, particularly, introduce the emergent mutual statistical interaction. 
In Sec. Ill we will formulate a general lattice field theory basing on the phase string effect and introduce a pair of 
unconventional order parameters. Using this general theory, we will present a general result for electric transport in 
Sec. |IV| In particular, we will prove a composition rule that expresses the electric conductivity in terms of holon and 
spinon conductivities. Armed with the pair of unconventional order parameters and the composition rule, we will 
explore the quantum phase diagram in Sec. |Vj (Throughout this paper, we shall focus on the underdoped regime.) 
Finally, we will conclude in Sec.jVl] Some of the technical details are regelated to Appendices [AfCl 



II. PHASE STRING EFFECT: A SHORT REVIEW 



For the self-contained purpose we shall start from a brief exposition of the exact phase string effect [H] for the t-J 
model (i, J > 0) 

H = t £ £ (4c,, + h.c.) + J J2 (Si ■ S t - ^) = H t + Hj (1) 

(id) a (id) 

on a square lattice, where h.c. stands for the hermitian conjugate. In this Section, in particular, we will discuss the 
physical picture of the mutual statistical interaction between holons and spinons, as well as its mathematical origin. 
In Eq. nil), c\ a (c ia ) is the (fermionic) electron creation (annihilation) operator on site i with the spin polarization 

a =4-5 T : n i = So- c lcr c ic 1S ^ ne occupation number operator, Si is the spin operator, and stands for the link 

between two nearest neighbors i and j. Note that throughout this work, we shall not distinguish notationally the 
difference between the operators and the numbers. The wave function is restricted on the Hilbert space such that 
no sites are occupied by more than one electron-the no-double-occupancy constraint. Eq. ([I]) together with the 
no-double-occupancy constraint serve as the exact starting point of the present work. 



A. Bosonization: phase-string representation 

In Eq. ([I]) the first (H t ) and the second (-Hj) term describe the charge hopping and the spin-flip process, respectively. 
The latter is reduced into the Heisenberg model in the undoped (half-filling) case where n, = 1 . For this case, long 
time ago Marshall showed that the matrix element of the superexchange Hamiltonian must satisfy a sign rule |15) . 
That is, there exists a complete set of spin bases {\4>}} such that given arbitrary spin configurations <p and <f>' , 

(<f>'\HM<0. (2) 

For doped antiferromagnets the Marshall sign rule |2| may be trivially realized within the Schwinger-boson (b\ a , bi a ), 

slave-fermion (//,/<) formalism, provided that the spin-hole bases {|</>)} are chosen appropriately. Under these bases 
the superexchange and hopping Hamiltonian are written as |I2j 

H * = -iEEW-A^- (3) 

<J<7' 

Ht = "* £ £ <rfUA*><° + h - c - ■ ( 4 ) 

(id) a 

Note that the sign of Hj is now negative. Here the "holon" creation operator fj and the "spinon" annihilation 
operator bi a commute with each other, and satisfy the no-double-occupancy constraint: 

///* + £ 4^ = 1- (5) 
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The sign structure in H t (arising from hole-spin exchange) results in peculiar properties. Technically, it prohibits a 
perturbative expansion of the resolvent: (E — Hj — H t )~ 1 even in t/J <C 1. Physically, it accounts for the fact [12l E] 
that the holon, upon hopping to some site occupied instantly by a f (i)-spin, acquires a sign +1 (—1). (As such, the 
Marshall sign rule is incompatible with the hopping process.) Therefore, two holon paths close to each other, though 
having (almost) identical amplitudes, may differ in the sign, because the parity of the number of hole-^-spin collisions 
of each paths may be different. This is the so-called phase string effect and possesses non-perturbative nature. 

Before further proceeding, it is necessary to take this non-perturbative phase-string effect into full account. Such a 
task was fulfilled in Ref. [H] by introducing the so-called phase string transformation: 

\ip) -> e ie e ze Oe~ ie , (6) 

which exactly "gauges away" the phase string arising from the holon hopping. In Eq. ^ the phase string operator 
is defined as 

6 = - E n'y e i3 n s u , % = Im \n( Zi — Zj) € [— k, +tt) , (7) 

where and n\ a are the holon and the spinon occupation number operator, respectively. Note that the two- 
dimensional lattice spans a complex plane and z is the position of the spinon (holon) in this complex plane. Moreover, 
the no-double-occupancy constraint implies that the summation excludes automatically the term with i = j . 

Let us now apply the phase-string transformation |6| to Eqs. (JsJ) and Q. Introduce the bosonic operator h\ , hi : 

h\ = fj cxp | -i E Q a n h x , hi= exp I i E B u n\ f t (8) 

which preserves the no-double-occupancy constraint: 

h\h i + J2blha = l. (9) 



Then, Hj t are transformed to 



-EE ^ b l b U^' A% . ( 1Q ) 

H t -> -t E ^ h t h J elaA ' h b )* b i° + h - c - . (11) 



(hi) a 



where the gauge fields A^f 1 given by 



A h = \ E (*« - 9 ») K ~ n h) mod27r . (i 2 ) 

A% = ^E (^-^)nfmod2^ (13) 

are angular variables. Eqs. (10)-(13) constitute the exact bosonization of the original t-J model ([I]) subject to the 
no-double-occupancy constraint. 

B. Mutual statistical interaction and compact gauge symmetry 



Consider now an arbitrary loop C on the lattice plane. Eqs. (12) and (13) give [T^] 



E A h = ^E ( 6 W- 6 W) mod2 ^< ( 14 ) 

(ij)ec ies c 

53 A% = tt E hjhi mod27r, (15) 

(i,j)ec ies c 
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where on the left hand side the sum runs over all the links on C , while on the right hand side the sum runs over 
all the sites inside C . Eqs. (14 1 and (15 1 show that the holon (spinon) carries a 7r-flux and couples to the motion 
of spinons (holons) via the gauge field vl^ (^tj)- We see that the phase string effect effectively "fractionalizes" an 
electron into two topological objects - the spinon and the holon - and the spin and charge degrees of freedom obeys 
a mutual statistical interaction. This kind of interaction was also found in different contexts [17H19| . 

It is of conceptual importance that the present theory differs from the earlier one [20] in the emergence of the local 
compact (instead of non-compact) U(l) ® U(l) gauge symmetry, where the two U(l) gauge degrees of freedom couple 
to the holon and the spinon, respectively. Indeed, Eqs. (10), (11), (14) and (15) are invariant under the following 



gauge transformation (All the quantities here may depend on the imaginary time r.): 



i,h _ 
h\ 



h\. 



(16) 



and are apparently invariant upon shifting A^ 1 by 2-7T . Importantly, due to the single- valued nature the following 
mapping 

from an arbitrary loop C s,h in spacetime lattice to the gauge group U(l) has the homotopy group Z. 



(17) 



////// 



////// 



WWW 



WWW 



FIG. 1: Dual lattice regularization: the lattice (red arrow) where spinons reside is separated from the lattice (blue circle) 
supporting the holon motion (presented by the path). 



III. LATTICE FIELD THEORY AND UNCONVENTIONAL ORDER PARAMETERS 



Our theory has so far been exact. In the remaining of this paper we will study the mean field approximation of 
Eqs. (flOb and (fill, read 



(ij) ° 

h h = i h y, < ' 1 i>:>>., ■ '>•<•• • 



(18) 
(19) 

It is 



The "minimal" Hamiltonian Hh + H s and Eqs. (12) and (13) constitute the so-called phase string model, 
important that this model preserves the compact U(l) <S> U(l) gauge symmetry. 

To further proceed we introduce the dual lattice regularization [20] , Basically, the spinon and the holon reside in 
their own square sublattice (Fig. [I]). The advantage of this ultraviolet regularization is as follows: as the constraints 
( 14 ) and ( 15 ) arc concerned, in the right-hand side of these two equations, the holon (spinon) does not appear on the 
loop C (We stress, however, that this is merely for the purpose of technical simplicity and violating this condition does 
not lead to changes of the results reported in this paper.). The spinon coordinate system is denoted as x = (ii, i2,io), 
where the first two components refer to the lattice point in the x-y plane and the last to the (discrete) imaginary 
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time. We set the lattice constant to unity. The holon coordinate system is denoted as x = (h, I2, Iq), with the origin 
as (—5,-5,0) in the holon coordinate system. The gauge field A x+ ~ x = A^(x) is defined on the spinon lattice, while 
A xx _. = A s ^(x) [as well as the external electromagnetic field A xx _^ = A^(x)] is defined on the holon lattice. Here, 
we denote the unit vector in the /1 direction as fi. Finally, wherever a lattice function f(i,I) is involved the relation 
(ii, = [i\ + 5 , i-2 + 5, *o) is implied for the arguments, and throughout this paper, the Greek indices fi 7 v, A run 

over both the spatial and the temporal components, while a,/3 only over the spatial component. 

Consider the partition function Z = Trexp(— (3%) . Here, j3 is the inverse temperature and the Hamiltonian is 

■H ee -t h J2 (e^h\hj. & + h.c.) - J s £ {e^b\ +&a b\_ a + h.c.) (20) 

la iacr 

under the dual lattice regularization. Then, it is a canonical procedure to divide /3 into TV pieces, and for each 
imaginary time slice we insert the unit resolution. Upon sending TV to infinity we expect to obtain an expres- 
sion of Z in terms of the path integral of the configuration of the gauge fields A s ^ h and the matter fields h, b 
(as well as their complex conjugates Jv,tf). All the fields are in turn imaginary time-dependent (which, however, 
is suppressed to simplify notations wherever no confusions arise). We now facilitate this general program by two steps. 

Step I. We absorb formally the constraints enforcing on the gauge fields by the compactness into the mea- 
sure, obtaining an expression for the partition function which is determined by an auxiliary Lagrangian C ; 

Step II. We work out the constraints explicitly, prompting £ to a new Lagrangian Jzf . 



A. Mutual Chern-Simons term 



We begin with Step I. To this end we note that for each imaginary time slice r = J/3/7V ,j — 1,2, ••■ ,7V, the 
topological constraints (IT4|) and (flBl imposed on each plaquette give 



& A . = e^d a A s p - 7T (bl f h t - fe] 4 6 u ) = . 
& A h = ^ af) d a A h ^ - nh\hj = , 



(21) 



where we define the (spacetime) lattice derivative as d fl f(x) = f{x +/t) — f(x) for the spinon lattice and as d^f(x) 
f(x) — f(x — p.) for the holon lattice, with fix) an arbitrary lattice function. Taking this into account we find 



D A [A s a ,A h a ]D[h\ h, b\b]8{,? A s)5(,? Ah )e 



-T, x c M \ 



A?=AZ=0 



(22) 



with the Lagrangian 



C M = h\(d -iA s + X h )hi + t h ^2(e iA "h\hi- a + c.c. 

a 

+ $>t(do - i*A h + A s )6 iff + J S J2 {^<b\ +&rT b\_ 



(23) 



Here c.c. stands for the complex conjugate, and the subscript "A" in the measure D A for that as the integral over 
the gauge fields are concerned, the constraints enforced by the compactness (which will be work out explicitly in the 
next subsection) are implied. The parameters X s ' h are to be determined below. The constraint S(JP A s) S(^ A h) may 
be released by introducing two real number valued fields Aq (x) , A$ (x) 6 K . As a result we obtain 



D A [A s a , A h a ] D A [A s ,A h ]D[h\h,b\b] exp 



C'c 



cs) 



where the mutual Chern-Simons term is 



cs 



- (A^d a A h + A h e°^d a A s ) 



(24) 



(25) 



The exponent in Eq. ( 24 ) is invariant merely under the imaginary time-independent transverse gauge transformation: 



4" 

A' 



A h a + d a t 



At 



he 



-i8 s 



(26) 



G 



with 9 s ' h some regular lattice functions. On the physical ground we expect that £' cs possesses a natural extension 
(221: 



£' cs -> £ns = -e^ x A s „d,A\ 



-cs 



(27) 



As such, the action: ~^2 x {Cm + £cs) is invariant under the full local U(l) ® £7(1) gauge transformation (16). In 
Appendix [A] we present an alternative derivation of this Lagrangian. 

Note that the elevation: £' cs — > Cqs is merely due to the lift of gauge fixing. This is most easily seen by reversing 
the procedure. Suppose that we choose the Coulomb gauge, the gauge fields are then composed of purely transverse 
components which may be explicitly written as [22] A s ± = {d y s , — d x 9 s , 4> s ) and A*\_ — {d y 9 h \ —d x 9 h , (f> h ) with cj) s ' h 
and 8 s ' h some functions in the spacetime lattice. (The first two components of A 8 ^ 1 are the spatial components while 
the last one the temporal component.) Inserting this expression into ^-CS we, indeed, end up with ^CS ■ One 
must caution that such an elevation may generally result in a normalization factor of Z accounting for the gauge 
degree of freedom. However, it does not lead to any physical results and we, therefore, shall ignore this factor. 

To proceed further, we divide the spinon lattice into two sublattices, say A and B. The spinon field in the sublatticc 
A is denoted in the same way as before while in the sublattice B is denoted as (6, b*). Furthermore, we introduce a 
four- component vector, (i g A), as follow: 



_ ( 1pi 



A = 



"4 



(28) 



Assuming the spinon fields are smooth over the (spatial) scale of the lattice constant, the Lagrangian Cm is 
expressed as 



re- 



-M 



(£> s + \ h - t h {2 + />)) h + , 



(29) 



M = 



-J s (2 + (D h )*) \ 

(D%)*+\ s -J s (2 + D h ) 

o -j s (2 + b h ) -(£#)* + x s o 

V -J s (2 + (/>)*) -D% + \ s J 



Here, Dg = d^ — iA^ and Dq = do — iAy are covariant time derivatives. D s = J2 a (^ a ^ ^a) 2 and D h = ^ a {d a — iAa) 2 
are (covariant) discrete Laplacians. (Notice that all the fields above are non-singular and, therefore, we may perform 
a hydrodynamic expansion. Alluding to this, the Laplacians result.) Finally, the partition function is written as 



Z = J D A [A s ,A h ]D[h\h,&,V] e -^^ £M+£ cs) 



(30) 



B. Lattice field theory 



So far, the construction of the field theory has been restricted on the topological trivial sector in the sense that 
the mapping (17) from arbitrarily given loop C s ' h onto 17(1) has zero winding number. Therefore, the field theory 



achieved in Step I does not carry any information regarding the compactness of the intrinsic gauge group. To fulfill 
Step II let us start from analyzing the mapping (17). It implies that in the spinon lattice, a loop C may pass through 
the area circulated by the loop C s once (Fig. [2| such that 



e» vX d v d x s {x) = 2wM s Y Jc'( x ')> M ' S e z 



(31) 



Here, the link field Jq(x) takes the value of +1 (—1) for the link (xx + fi) ((x + fix)) on C and zero otherwise. Under 
the gauge transformation, the action *^2 x (£-m + £cs) is transformed to 



^(£ M + C CS ) -> + C CS ) + ^ uX d,0 s d u A h x = ^(£ M + C CS ) - 1 Y ^ A h x d^9 s , 



(32) 



where in deriving the second equality we used the integral by parts. In order for the lattice field theory to be unaffected 
we demand that the second term of Eq. p2l to be multiple of 2it . Taking Eq. (31 ) into account we obtain 



Y J£(x)4(x) = 0mod7r. 



(33) 
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Using the Stokes' theorem, we obtain Y^ x (dfiA^ — d u A^)AS^[Sc] — Omod n with ASg u [Sc>] the unit surface element 
at x . (Note that here no summation is implied for /x, v .) Since the loop C and thereby the surface Sc are rather 
general, Eq. (33) then implies 



d^A^ - d u A'^ = mod n . 



(34) 



That is, the Maxwell tensor of A is locally quantized, with a quanta it . Likewise, we have 



d^Al - d v A s = mod 7r . 



(35) 



c 



c 



FIG. 2: The loop C passes through the area circulated by another loop C s 



Eqs. ( 34 1 and ( 35 1 are the exact constraints absorbed in the measure D A [A s ,A h ] of Eq. 30}. To release such 
constraints we further introduce two dual integer fields JV* L and . In doing so we eventually obtain the path 
integral representation of the partition function 2, , which is 



Z 



f D[A s ,A h ]D[h\h,¥,^}expl~^2[C M 



£ 



CS 



2te^ x (A h ^ x s + AfrJ?)] 



(36) 



C. Symmetries 



We turn now to examine various symmetries of the action in Eq. (36). First of all, it is apparently invariant 
under the local U(l) (g> i7(l) gauge transformation (16). Then, we absorb the integer fields ,jV s ' h into the mutual 



Chern-Simons term. In doing so we reduce the Lagrangian in Eq. (36) into 



C M + - e^ x (A* - 2^) d v (A h x - 2n^ x h ) , 



(37) 



where we ignore the irrelevant term multiple of 2n. From this we see that the field theory is compact in the gauge 
degrees of freedom, because a multiple 2ir shift of the gauge fields is absorbed by the integer fields. Such property 
is intrinsic to the discrete lattice geometry but not to the dual lattice regularization. It is also easy to see that the 
Lagrangian possesses both parity and time-reversal symmetry |20) . Under parity transformation jV h transforms as 
an axial vector while jV s as a polar vector. Under time-reversal operation ,yV s ' h transforms as 



A, 8 



A/a A/s 



"'a i 



(38) 



We further show that the Lagrangian possesses the spin rotation symmetry. Indeed, transforming as 



e -^*"/ 2 



U{- 



-4>,-x) 











(39) 



we find that the Lagrangian is invariant. Note that the 2x2 matrices above _are defined on the sector introduced 

10' 



by the doublet due to the sublattice structure [cf. the definition of 4^ in Eq. (28)], while U and a 3 = 



-1 



arc- 



defined on the spin sector [cf. the definitions of ipi , 4*i m Eq. ( 28 )] . The spin rotation is generated by the Euler angles 
0, <p and x, i.e., 

rrm j. ^ ( cosfe^/ 2 e^/2 8 i n 4 e -**/2 e <x/2 \ 

"(CM = [ _ sin | e#/2e -i X /2 cos | e -^/2 e -ix/2 J • (4°) 



Finally, the single- valueness of transformed "P; is guaranteed by the relation d^ — 2A^ and Eq. (34 1. The present 
theory shows that the spin rotation symmetry is protected against high-energy ferromagnetic fluctuations, while the 
earlier theory |20j proves the existence of this symmetry only for the low-energy sector. 

To proceed further, we wish to soften the hard-core boson condition, i.e., Eq. ([9|. In doing so, all the above 
symmetries must be respected. This goal can be achieved by modifying the Lagrangian to be 

S£ = C h + C s + i e^ A {A% - d v (A h x - 2ttJK x h ) , (41) 

where the first two terms are the spinon and holon Lagrangian, respectively, read 

C h = /;t (pi + A* - t h (2 + !>)) h+ f {tftif , 

with the last term in s describing the on-site repulsion (the coefficients ui 3> tf ll u 2 3> J s ), while the last term in 



_Sf [cf. (41 1] is the compact mutual Chern-Simons term. Thus, the partition function is promoted to 



Z= D[A s ,A h ]D[h\h,^,^}e~ s (43) 

with S = ^£ . Finally, X s ' h are determined by the minimum of the effective action obtained by integrating out 



(sum up) all the fields. Eqs. (41 )-(43| complete the construction of our compact mutual Chern-Simons theory. Notice 
that here A^ h are compact degrees of freedom, i.e., A s J l £ [—it, +tt], while Aq H e M. not. 

It is important to remark that the integer field jV£ (^a) captures the singular part of the phase fluctuations of 
the spinon (holon) field. As we shall see shortly later, in the presence of the spinon (holon) superfluid ,jY£ {-V^a) is 
indeed the integer field introduced in the well-known Villain's approximation [5T]. However, the integer field jVq' 
have completely different physical meaning: (i) Summing up JV^ leads to the quantization of e 0a Pd a Ap with a unit 
7r . By integrating out Aq , we then find Wh G Z . Thus, in the limit u\ — > oo , the integral over the holon field is 
dominated by the holon configuration where |/i(a;)| 2 takes the value of or 1 . (ii) Similarly, summing up jYq leads 
to the quantization of e 0a ^d a Ap . By integrating out Aq , we then find ^/tj]^ g Z/{0} , where S = diag(cT3, a?,) T is 
defined in the sector introduced by the sublattice structure. Thus, in the limit u 2 — > oo , the integral over the spinon 
field is dominated by the spinon configuration where for given x, either |&^(x)| 2 or |6^(a;)| 2 takes the value of unity. 
Together with the dual lattice regularization, (i) and (ii) realize the no-double-occupancy constraint (p)J). 



D. Unconventional order parameters 



The Lagrangian (41 ) shows that the two gauge fields, A 3 ' , are dual but the matter fields \E' and h not. The duality 
of the gauge fields suggests the introduction of a pair of Wilson loops, defined as 

W' h [C] = Z- 1 I ' D[A\A h ]D[h\h^\^]e- s+l ^* A » hj c , (44) 

which depends only on two parameters namely the temperature and the doping. Here, C is a spacetime rectangle with 
length T (R) in the imaginary time (spatial) direction. 

Physically, the Wilson loop Vy s [C] (W [C]) probes the interaction of a pair of test holons (spinons) at a distance 
R, V h (R) (V s (R)), via 

V h (R) = - limbec f In W s [C] , , , 

V s (R) = - lim T ->.oo i In W h [C] . [ ' 

Furthermore, the analysis in the remaining of this paper suggests that this pair of Wilson loops suffices to characterize 
the phase diagram of the t-J model. Therefore, it plays the role of "order parameter" and informs the non-Landau- 
Ginzburg- Wilson nature of phase transitions involved. We remark that the Wilson loops introduced here differ 
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crucially from a conventional one denned on a pure gauge field theory. In fact, it is the coupling between the matter 
and the gauge degrees of freedom that leads this pair of Wilson loops (potentially) to display very rich behavior. 
Physically, the existence of the Wilson loops as order parameters is the reminiscence of the strategy adopted in 
quantum chromodynamics, where the Wilson loop serves as a canonical order parameter to distinguish the cold and 
plasma phases [2"5j . 

IV. COMPOSITION RULE FOR ELECTRIC TRANSPORT 



In this Section, we will present a microscopic theory of electric transport basing on the lattice field theory (41 1. 
Specifically, we will derive a so-called composition rule that expresses the physical electric conductivity tensor, a e = 
{cr" in terms of the holon conductivity tensor, a h = {a" 13 }, and the spinon conductivity tensor, a s — {a" ^}. We 
will begin with a phenomenological discussion on this rule and then proceed to a microscopic justification. 

A. Phenomenological discussions 

The minimization of the effective action with respect to gauge fields leads to the following equations of motion: 

4; = % + ^" X d v (A' x - 2*JT X ') = 0. [b) 
where the spin and charge currents are defined as j s / h = innr/i ■ From these equations we obtain 

■h 1 ,0afl T?h -s _ 1 Ja^ ps /A7\ 

Ja ~ - £ Ja ~ -£ (47) 

7T 7T 

where we have introduced the (macroscopic) internal "electric" fields in the imaginary time representation, 

E s a ' h ee -ie^d^ (A s u - h - 2T:,yV l f' h ) . (48) 

On the other hand, in the presence of an external electric field E e a which couples merely to the holon degree of freedom, 
the linear response assumes 

o h M = ofbWpiz) + m% fail) = *f(q)E h e(Q)> ( 49 ) 

where we have passed to the Fourier representation with q = (q x , q y , iuj n ) the Fourier indices and uj n the bosonic 



Matsubara frequency. Combining Eqs. (47 1 and (49), we find that the electric conductivity tensor, defined by 

& = o?E%, (50) 

obeys the following composition rule: 

a e l =a h 1 - n 2 ea s e , (51) 

where e = {e Q ^} is the antisymmetric matrix, with e xy — —e yx = 1 . Note that here and after, to make formula 
compact we shall omit the argument q wherever confusions may arise. 

B. Microscopic justifications 



Now we turn to present a perturbative proof of the composition rule (51 ). On this purpose we note that in each 
phase one may in principle integrate out all the matter fields, arriving at an effective action of pure gauge fields. 
Because the holon (spinon) field couples only to A s (A h ), this gauge field action has the general structure, 



S^A 8 + A e ] + S 2 [A h ] + - V e^ x (A s - 2-KjY*)d v (A\ - 2-KjV x h ) (52) 
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and is gauge invariant, where the first two terms are the kinetic part. We apply the gauge fixing condition Aq = Aq = 
Aq = and expand the effective action around its saddle point in terms of the fluctuating gauge fields A^ h . Keeping 
the expansion up to the quadratic order, we obtain the fluctuating action, 



{[A' a (x) + A e a (x)]uf(x - x')[A^(x') + Ap (x 1 )] + A*(x)ILf(x - x')A h p {x')) - % - £ ^A s a d A h p , (53) 



where IT^ iS are the polarizations of A s and A , respectively. Integrating out A s ' h eventually reduces the partition 
function to Z oc e~ Scff ^ A 1, where the prefactor is independent of A e and iS e ff[A e ] is the effective action of external 
gauge fields expanded to the quadratic order of A e : 



(54) 



Here, II e = {I!"' 3 } is the polarization of A e , expressed in terms of ±1^ = {n"^} via 



n„ 



2 

7T = 

— eIT s e. 



(55) 



Noting a^f 



e, h, s), we obtained from Eq. (55 1 the composition rule (51 ). 
For later convenience, here we consider a simplification of the combination rule (51) by ignoring the off-diagonal 



components of the conductivity tensor (namely the crossing transport). At zero temperature, the static conductivity 
may be obtained by taking the limit, q a — > 0, first and then uj — > 0. For an isotropic system the conductivity tensor 
is reduced to <j% 



' x 



&x5aBi an d the composition rule (51 1 is reduced to 



7T 2 CT S 



(56) 



We further provide a qualitative explanation of this rule. The macroscopic electric current (density) j\ is fully carried 
by holons and driven by both the external electric field E a and "electric field" E^ induced by spinons. The latter 
finds its origin analogous to that of Ohmic dissipation in type-II superconductors: each spinon mimics a "magnetic 
vortex" suspending in holon fluids and, upon moving, generates an elect ric field antiparallel to j^, i.e., E s a = — 7r 2 cr s j'^. 



In combination with the Ohm's law, i.e., a h = — n 2 c7 s j^ + E a , Eq. (|56[) then follows. 



V. QUANTUM PHASE DIAGRAM 



In this Section we shall consider an application of the general theory developed in Sec. Ill and IV Specifically, we 
will calculate the Wilson loops, W^'^C], at zero temperature where such a pair of non-LGW order parameters depends 
merely on the doping. We will show that the order parameter FF S [C] (W [C]) displays non-analyticity where a holon 
(spinon) deconfincmcnt transition occurs. This pair of non-LGW order parameters indicates that the quantum phase 
diagram (see Fig. 3) is composed of three phases, namely the antiferromagnetic phase (adjacent to zero doping), the 
superconducting phase (far away from zero doping) and a novel phase - the Bose-insulating phase (for intermediate 
doping) . 



A. Superconducting phase 

We first study the regime far away from zero doping, where the superconducting phase is formed. In this phase, the 
spinons are confined while holons deconfincd. Such peculiar confinement properties of spinons (holons) are intrinsic 
to the onset of holon superfluids. Upon decreasing the doping, the spinons undergo a deconfincmcnt transition. 
Correspondingly, the system loses its superconducting long-ranged order. 



1. Logarithmic spinon confinement and holon deconfinement 

In this regime, the holon Lagrangian Ch is minimized upon the formation of superfluid, justifying the supercon- 
ducting phase (with a density nh — (2th — A' l )/ui). In the presence of a holon superfluid, a uniform "magnetic field" 



11 



IK**** 




<3 



<3sfc 




antiferromagnet 
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FIG. 3: (Color online) Quantum phase diagram of underdoped Mott insulators. The ball (in blue) and the arrow (in red) 
stand for the holon and spinon, respectively. The vortex, in red (blue), surrounding a holon (spinon) arises from the spinon 
(holon) condensate. The wavy line stands for the confinement. 



is formed and minimally couples to the motion of spinons [cf. Eq. (42)]. Consequently, a Landau- type gap is opened 



in the spinon spectrum and spin excitations are suppressed. For sufficiently large doping, these spin excitations may 
be ignored and the functional integral over the spin degree of freedom *f> is frozen at state composing of the so-called 
resonance valence bond (RVB) pairs. As such, the functional integral over iff accounts only for an irrelevant overall 
factor, and the Wilson loop, W-^C] , is simplified to 



W h [C] oc / D \- A *' A h ]D[tf, h] expl-J2 



C h + V" A (A; - 2irJ£) d v (A h x - 2nJK x h ) - U*Jg 



(57) 



To proceed further, we separate the spatial components of the gauge field, A\ — 2-K^y^ , into two parts. The first, 
denoted as A' ( ' t , is the background component which satisfies 



<? afi d n A\ = 



a Ap = Trn h 



(58) 



and therefore, is imaginary time-dependent. (To be consistent with this condition it is necessary that .Aq s must be 
uniform in the lattice plane.) The second is the fluctuation component. Then, we factorize the holon field as 



= IUp- ifl iW 



h(x) = \h\e 



(59) 



substitute it into Eq. (57) and make the change of variable: A 8 a — > A s a — d a 9\ . Integrating out the matter field h with 
the help of saddle point approximation and summing up the integer field JV^ 1 , we obtain 



W h [C] oc / D[A',A h ]exp\-Y, 



2*JT') d„A h x + -e 0ap A s d a A p 

7T 



exp < i 



x exp 



{-^X {A ° +2nhthUiA%2) } 



(60) 



Here, the integral variable A\ stands for the fluctuation component of the corresponding gauge field. In deriving this 
equation we have used the facts that summing up the fluctuating flux gives a vanishing result, and that ,jV$ is uniform 
in the lattice plane. Notice that an irrelevant factor e % ^^ A c J c has been omitted. Importantly, the first term in the 
first exponent above suggests that in the presence of holon superfluids, Jf£ is identical to the integer fields introduced 
by the Villain's approximation [3T]. Notice that d\ 1 d^9\ 6 [— n, +tt] is non-singular, and the singular component of 
the phase fluctuations is very captured by (the spatial components of) the integer field. Integrating out the A s fields 
gives 



W h [C] (x ]T J D[A h ] exp J - Y, fe (^X) 2 + (^"d^Ai) 2 - iA\ (j£ + 2^ a d v ^) 



(61) 



From this we see that the dynamics of the fluctuations of the gauge field A h displays emergent Lorentz symmetry 
with a "speed of light" ci = \j2ny l ty l u\ . By introducing the following rescaling: 



(d / Cl ,d x ,d y )^d^ {A%/ Cl ,A h x ,Al)^A% (ciJq, Jq , Jq) J[ 



C ' 



(62) 
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we rewrite Eq. (61) as 



W h [C] oc J2 I 


f D[A h ]exp< 












f D[A h ]exp< 


[-E 









1 

4ef 

— f-F fc 

4e 2 v M« 



z/l£ j£ - 2i^e OQ/3 d a ^ a s 



(63) 



Here, F^ v = d^A^ — d v A^ is the Maxwell tensor and ef = n 2 /ui is the squared "bare charge". In obtaining the second 
line we have used the fact, d^jV^ = , at zero temperature. 

If no phase vortices are present, i.e., 2ire Qal3 d a tA / p s — 0, Eq. (63 1 is simplified to 



W h [C] oc J D[A h ]exp- 



E 



4ef 



2 \ A iiv) 



As we will show in Appendix [B] it can be further reduced to (Here, we restore the original unit.) 

R 



W h [C] oc e - Tys{R \ V s (R) 



(e lCl )\ R 

— 7. In. — 

2tt 



irrihth In - 



(64) 



(65) 



for T ~ > oo, where R{ is the ultraviolet cutoff. Notice that the coefficient of the potential, i.e., irrihth , does not 
depend on the strength of the on-site repulsive interaction. Eq. ( 65 ) suggests that an external dipole undergoes 
logarithmic confinement and has important consequences: (i) The presence of a pair of free phase vortices of the 
holon superfluid, carrying opposite vorticity, is energetically unfavorable which is consistent with the simplification 
above; (ii) In the presence of spin excitations (right, Fig. [3]), this dipole may mimic a pair of spinons with opposite or 
identical polarizations. In the latter case, a phase vortex with a vorticity of ±2-7r may be excited from the background 
and bound to a spinon of flux =p7r, reversing the sign of the bare topological charge of the spinon accordingly, i.e., 
=p7r — > ±7r . That is, a pair of spin excitations must be logarithmically confined. 

We turn now to calculate W s [C] . To this end we may ignore the fluctuation component of the gauge field A\ . 
Following the same procedures of deriving Eq. (60) from Eq. (57), we find (Here, we do not rescale the gauge field 
and d^ .) 



W S [C] 



D[A s ,A h ]ex V \-J2 



2u x 



(A s 2 + n h t hUl A s a 2 ) + 



f 0a/3 



(A s a -2nJK a s )d^ 



(66) 



where we have used the fact that 9\(x) is non-singular giving d a dp9i = dpd a 6\ and d^OiJ^ = . The functional 
integrals over the temporal and spatial components of A s are factorizable, i.e., 



W S [C] oc J D[Aq] exp j - £ - *W 



exp 



D[A' a) A%] CX P - 



E 



n h t h A s a 2 + ~e a0 ° (A s a - 2ir,yV a s ) A h Q 



fl 



y/nhth 



+ h 



R 



V n hth 



(67) 



with the scaling function fx (/ 2 ) independent of R (T). This is none but a (holon) deconfinement law - insensitive to 
the details of the scaling functions fi t 2 - as expected by holon superfluids. Shortly we will see that this Wilson loop 
is not critical at the spinon deconfinement critical point. 



2. Non-BCS nature of superconductivity 

Let us subject the system to a small magnetic field with a total flux <I> e . Consider an area A on the spatial plane 
that is enclosed by a sufficiently large loop. Since in the ground state, spinons are confined into RVB pairs and do 
not contribute a net flux, integrating out the Aq field leads to 

Y * f S* = $ e mod 2vr . (68) 
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As we explained above, the left-hand side is quantized with a quanta n . This implies in turn that Eq. ( 68 1 is none 
but the quantization condition of the external magnetic flux, i.e., $ e = n(hc/2e), n £ Z in the full unit |16j. However, 
this quantization is intrinsic to the topological spin excitations and differs in the nature from its counterpart in BCS 
superconductivity. In fact, in the absence of the external magnetic field, i.e., $ e — , a single spin carrying a flux 
of ±7r cannot be excited otherwise Eq. (68) is violated. Instead, spins are excited in pairs, constituting a spin-0 (or 
spin-1) excitation. However, in the presence of the external magnetic field, i.e., <I> e ^ , a single spin can be excited 
provided it is nucleated at the magnetic vortex core. Moreover, as a result of the spin rotation symmetry, both 
polarizations are possible [2~3] . 

We then examine the static electric conductivity. As mentioned above, in the absence of the external magnetic field, 
spinons are excited in pairs. If two spinons have identical polarization, they cannot mobile because they are always 
bound to a (local) phase vortex of holon supcrfluids, as shown in Sec. V A 1 If two spinons have opposite polarization, 



as shown in Sec. |V Al) they undergo logarithmic confinement. Therefore, neither of these spinon pairs supports spinon 
transport, i.e., a s = 0. On the other hand, since holons undergo Bose condensation, the holon conductivity is infinite, 
i.e., ah — co. From the composition rule (56) we then obtain a e = oo. It is important that according to composition 
rule, the establishment of superconductivity relies crucially on a vanishing spinon conductivity a s which arises from the 
spinon confinement. In other words, it suggests that the disappearance of superconductivity has to be associated with 
the spinon deconfinement, where a s no longer vanishes - a fact reflecting the non-BCS nature of superconductivity 
that we will see in Sec. IV CI 



B. Antiferromagnetic phase 

Now we study the regime adjacent to zero doping, where the antiferromagnetic phase is formed. In this phase, the 
holons are confined while spinons deconfined. Such peculiar confinement properties of spinons (holons) are intrinsic 
to the onset of spinon superfluids. Upon increasing the doping, the holons undergo a deconfinement transition. 
Correspondingly, the system loses its antiferromagnetic long-ranged order. The analysis below is largely parallel to 
those of Sec. |V A| Therefore, we shall only sketch the main steps. 



1. Logarithmic holon confinement and spinon deconfinement 



Consider the case where holons are sufficiently dilute and the holon Lagrangian is thereby ignored. As a result, 



(69) 



The antiferromagnetism is justified by the existence of homogeneous saddle points, denoted as ,i- e -j 

= . (70) 



SS 



In AppendixK^lwe will show that (i) the saddle point has the structure as \&o = d/^? 6 *^' y/™l elSi '■> \f^i el9i '> -v/% e * 0t ) T 
with n a and ^homogeneous in spacetime and n^+n^ = (2J S — X s )/(2u2) ; and (ii) the total spin polarization vanishes, 



as a manifestation of antiferromagnetism. 

Next we wish to integrate out the spinon field. To this end we factorize the field as 



(71) 



Sn 1 (x)+Sn 2 {x) i6 9 ,(x)/2 \ 



^l n ^ | &n 3 (x)-5m(x) c -i9 9 (x)/2 



6n 3 (x)+&n l (x) i8 2 (x)/2 



(72) 



^ ^J nt + Sn 1 {x)-8n 2 (x) e -i9 2 (x)/2 ^ 
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and insert it into Eq. (69 1. The phase field 6*2(0;) generates the Goldstone mode. Furthermore, because of Eq. (71) 
the background component of the gauge field A s a — 2-xjV£ vanishes. With the help of the saddle point approximation 
we obtain from Eq. ( 69 ) 



W S [C] oc ]T [ D[A s ,A h ]D[6n,e 2 }exp\-J2 -e avX {A h a ~ 2^^) d u A s x + -e QaP A h Q d a Al 



x exp 



+ 2J s (n t + n_i)(d 



h\21 



expi^ 



J* 



{5n\ + 8n\ — 5n\ — Snfj 



y {8ni + Sn 3 ) 2 



exp < 1 



(73) 



Again O2 (x) is non-singular and instead, the integer field JV£ characterizes the singular part of the phase field namely 
the phase vortex of the spinon superfluid. Integrating out Sn and 9 gives 



W*[C] ex ]T / D[A\ A s ] exp <j - £ ^ (Af + W(n t + n^A™) + % - (A\ - 2^^) d v A\ - za* j£ 



(74) 



This result is identical to Eq. (60 1 upon making the replacement: u 1 1 — > J ' s / '4, 2rt/ l t/ ! .Mi — > 16(n^ + n^) for the 
parameters and h — > s, s — > h for the superscripts. Translating Eq. (65) into the present context, we find 



W S [C] oc e - Tvh W, V h (R) = 27r(n t + n 4 ) J s In 



R 

tin 



(75) 



for T —> 00 . Here, i?2 is the ultraviolet cutoff. Notice that here the Lorentz symmetry emerges but with a different 
"speed of light" . The squared "bare charge" is also different from that of the superconducting phase. 

Eq. ( 75 ) also suggests that an external dipole undergoes logarithmic confinement and has important consequences 
as follows, (i) The presence of a pair of free phase vortices of the spinon superfluid, carrying opposite vorticity, is 
energetically unfavorable, (ii) In the presence of holon excitations, a phase vortex with a vorticity of —2tt may be 
excited from the background and bound to a holon of flux +7r, reversing the sign of the bare topological charge of the 
holon accordingly (the so-called "anti- holon" ) , i.e., tt — > —n. Eq. (75) then implies that such a holon-anti-holon pair 
is bound together via the logarithmic confinement (left, Fig. [3]). In other words, two holons are bound to a phase 
vortex of the spinon superfluid with a vorticity of — 2n . 

We turn now to calculate W h [C] . Similar to Eq. ( 66 1 , we find 



W h [C]<x 
which gives 



D[A\ A s ] exp { J2 T ( A ° 2 + 16(n ^ + ^ )A " 2 ) + I ^ t A " ~ d ^ - lA * J ' 



W h [C] cx exp 



h 



c 2 T 



R 



,(76) 



(77) 



It suggests spinon deconfinement as expected by spinon superfluids, insensitive to the details of the scaling functions 
fx 2 ■ As we will see below, this Wilson loop in non-critical at the holon deconfinement critical point. 



2. Unconventional antiferromagnetism 



Similar to the discussions in Sec. VA2 the antiferromagnetism (<5 ^ 0) here is unconventional. Let us subject the 

Since in the ground state, holons and anti-holons are 



system to a spin twist generated by an external gauge field A' a 
confined in pairs and do not contribute a net flux, integrating out the Aq field leads to 



\ " Oa/3 



e Oati d a A' mod27r 



(78) 



dual to Eq. ( 68 1 . Since the left-hand side is quantized (with the same quanta 7r) , this implies a dual quantization 
condition: the external flux generating the spin twist may penetrate into the antiferromagnet only if it takes the value 
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of rar, neZ. This quantization is intrinsic to the topological holon excitations. In fact, in the absence of the external 
spin twist, a single holon cannot appear in the excitation spectrum otherwise Eq. ( |78[ ) is violated. Instead, the holon 
and the anti-holon are excited in pairs, constituting a charge-2 bosonic excitation. However, in the presence of the 
spin twist, a single (anti-)holon can be excited provided it is nucleated at the center of the spin twist. 

For the static electric transport, we note that in the absence of external spin twist, the holon and the anti-holon 
are excited in pairs. According to Sec. |VB 1| such pair is bound to a phase vortex of spinon superffuids. Since the 
latter is localized in space, (As such, the spontaneous translational symmetry breaking appears.) — . From the 
composition rule (56 1 we then find that the antiferromagnetic phase is insulating, i.e., a e = . It is important that 
such a property of electric transport is intrinsic to the holon confinement and, therefore, persists to some finite doping 
- the holon deconfinement critical point. This is in sharp contrast to previous theories [6] where the superconducting 
phase was found to be pushed all the way down to zero doping. Finally, it should be noted that without the integer 
field jV^ describing the spinon phase vortex, such an antiferromagnetic insulating phase cannot be established. 



Bose insulating phase 



The qualitatively different behaviors of W [C] described by Eqs. ( 65 1 and ( 77 ) suggest a critical doping 82 , at which 
W [C] - as a function of S - is non-analytic. This is the spinon deconfinement quantum critical point. Likewise, the 



expressions of Eqs. (67) and (75) for M^ S [C] suggest another critical doping <5 1; at which iy s [C] is non-analytic. This 
is the holon deconfinement quantum critical point. The mechanisms underlying these two quantum critical points - 
the disappearance of the superconducting or the antiferromagnetic long-ranged order, are independent. Thus, d± =/= 62 
generally. Since the confinement of holons (spinons) is a consequence of spinon (holon) condensation, the case of 
Si > 62 is ruled out. That is, spinons and holons cannot be confined simultaneously. Instead, we have 5% < 62 
generally. As such, there is an intermediate phase separating the antiferromagnetic and the superconducting phases. 

Indeed, this is a novel phase where both matter fields undergo condensation characterized by the field ^0 {ho)- It 
should be contrasted to the antiferromagnetic (superconducting) phase where only the spinon (holon) field 'J/ (h) is 
condensed. The fields vPq > ^0 are determined by the following self-consistent equations: 



= 0. 



ho, d o h o =0 



ss 

<5vfT 



= 0. 



*o,do*o=0 



e 0af3 d a (A' p - 2ttJK s ) = tt* |,Etf • 



(79) 

(80) 
(81) 



In general, the field ^0 {ha) has an amplitude inhomogeneous in space. In the Wilson loops (44) the functional integral 
over A s a , jV£ {Jvt , <J^a) i s dominated by the small fluctuations (denoted by a^ h ) near the 
satisfying Eqs. (80) and (81) 

1 j 1 AS.h s,h 

compact we denote A ■• ■■ 



Dackground configurations 
Taking Eqs. (79)-(81) into account, we rewrite Eq. (44) as (To make the formula 



as a, 



W s ' h [C] oc J D{a s ,a h ]D{h\h,¥,*}expl-Y, [c h \ a s- ~ ia^ h - h%) 



exp J -]T [C s \ <=0 - «4(*t£* - *t E$0 ; 



exp 



E 



^ vX n s rl n h — in s ' h 



(82) 



The first (second) exponent involves merely a s (a h ). It is important that integrating out the holon (spinon) field leads 
to an effective action of a s (a h ) which is massive. (This can be readily seen by observing that a spacetime- independent 
a s (a h ) cannot be absorbed into the functional integral over h (vt) and as such, the effective action of a homogeneous 
a s (a h ) must not vanish.) Therefore, the residual mutual statistical interaction between gauge field fluctuations, of 
order 0{du) , is negligible, and Eq. (82) is simplified to 



W s < h {C] oc J D[a s ,a h ]D[rf,h,¥,*}expl.-J2 [U\a^o - ia^ h - h%) 



x exp J - [£.|«s=o - ^0 (* f £* - ^S* ) 



exp 



(83) 
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Integrating out the matter fields then gives 

W s ' h [C] cx J D[a s ,a h ] exp j - ^ 



±- (a? + ( Cl <) 2 ) + f (af + (c 2 a h a f) - iafj£ 



(84) 



From this equation we see that the functional integrals over the temporal and spatial components of a s (a h ) can be 
factorized. As a result, we obtain 



w ' ,[ciot = p H A (J§ 

W h [C] oc exp 



thj \V n hth 



c 2 T 



R 



(85) 



where the new scaling functions /i j are independent of i? (T). Insensitive to the details of these scaling functions, 
the Wilson loops ( [85] ) indicate that an external pair composed of holon and anti-holon or spinons with opposite 
polarizations undergoes deconfinement (middle, Fig. [3]). Comparing these two expressions with Eqs. (67) and (77 1, 
we find that iy s [C] (W-^C]) is indeed non-critical at the spinon (holon) deconfinement critical point. 

The spinon and the holon deconfinement have far reaching consequences. First, both the long-ranged antiferromag- 
netic and superconducting order no longer exist in this phase. Second, recent studies by one of us |26j have shown 
that due to the spinon and the holon condensation, the phase of the electron operator is disordered, indicating the 
existence of a disordered gapless fermionic mode. As a result, this phase is compressible. (The present Bose insulator 
may be considered as a new example of the zero-temperature compressible quantum matters proposed recently |27j.) 
In other words, the holon (charge) density may be continuously tuned from Si to S 2 . Finally, alluding to a s ] % = oo 
arising from spinon (holon) condensation, we find that this phase is also insulating, a e — , from the composition 
rule (56). 



VI. CONCLUSIONS 



The t-J model (on a bipartite lattice) displays a non-perturbative sign structure, the so-called phase string effect. 
With this effect being taken into full account an electron is necessarily fractionalized into two bosonic constituents, the 
holon (the charge degree of freedom) and the spinon (the spin degree of freedom). Each constituent is a topological 
object and carries a tt flux. The latter mediates a compact U(l) gauge field, A s ^ (A'*), minimally couples to the 
motion of holons (spinons)-the so-called mutual statistical interaction. In this work, the exact phase string effect 
is refined in terms of the lattice field theory. Based on this field theory, a pair of unconventional order parameters 
namely the Wilson loops W s, ' l [C] is introduced. These two unconventional order parameters describe the holon 
(spinon) confinement-deconfinement property and suffice to characterize the phase diagram of the t-J model in the 
undcrdoped regime. We further establish a general composition rule for the electric transport, which expresses the 
electric conductivity in terms of the spinon and the holon conductivities. 

The lattice field theory and the general composition rule are applied to study the quantum phase diagram. In 
the antiferromagnetic regime, where the doping is sufficiently small, spinons undergo deconfinement while holons 
are confined, leading to an (electrically) insulating phase with long range antiferromagnetic order. Whereas for 
sufficiently large doping, holons undergo deconfinement while spinons are confined, leading to a superconducting 
phase. (Further analysis of such a superconducting phase unveils that this is a d-wave superconductor and possesses 
fermionic Bogoliubov qusiparticles [26 .) We find that W^fC] (IF^C]) displays non-analyticity at some doping 6i 
(62 with Si < 5 2 ), where the system undergoes a holon (spinon) deconfinement quantum phase transition. Most 
strikingly, we find that between the antiferromagnetic and the superconducting phase, there is a novel phase, the Bose 
insulating phase. In this phase, despite of spinon and holon condensation, no long range orders occur and the system 
is also electrically insulating. These results inform profoundly non-Landau-Ginzburg- Wilson nature of quantum phase 
transitions in the t-J model. We remark that different from the earlier field theoretic formulation |20j . the present 
lattice field theory is compact, and the compactness of the emergent U(l) gauge fields is essential to the formation 
of the antiferromagnetic and the Bose insulating phase. Finally, we should emphasize that the present theory is not 
limited to the quantum case. Applications to the finite temperature case are of fundamental importance and of great 
interests, which we leave for future studies. 

Note added. After this paper was submitted, we became aware of the works by Tesanovic and co-workers [28] - 
These authors also found that the duality between the superconducting and the antiferromagnetic phase dooms to 
lead to an intermediate insulating phase. Although this scenario is similar to the quantum phase diagram presented 
here, the intermediate phase reported in these papers has very different nature. 
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Appendix A: Alternative derivation of mutual Chern-Simons term 



In Sec. Ill A in deriving the mutual Chern-Simons term the gauge fields A s a ( A 1 ^) are essentially considered as 
"external" parameters as the holon (spinon) motion concerned [cf. Eqs. (18) and p9|]. These external parameters 
are then subject to the self-consistent constraints (21). In this Appendix, we follow Ref. l20l to give an alternative 
derivation of the mutual Chern-Simons term. 

The starting point is the observation that A s ^ h constitute the canonically conjugated pair. Indeed, with the help of 
the constraint ,^A h — [cf. Eq. (21)], the holon current conservation law may be rewritten as 



Jafl 



SH 



(Al) 



Exchanging the order of the derivatives of the first term, we obtain 

1 



(d x d A h -dyd A\ 



5H 



SH 

~dA°, 



(A2) 



which gives 



1 d A h - 5U 



sn 



-doA^- 

7T OA 



(A3) 



From the spinon current conservation law we may obtain two similar equations. Together with Eq. (A3) they justify 
the canonically conjugated relation: 



r\ h — A s 

ll x — A y ! 



Hy — — A S 



(A4) 



Passing to the path integral representation, instead of Eq. (24), we obtain Eq. (30) directly. 



Appendix B: Derivation of logarithmic confinement 



In this Appendix we prove Eq. ( 65 ) . Throughout this Appendix we set c\ = 1 , and without the loss of generality, 
we have set the four corners of the timelike rectangle to be (0, 0, 0), (i?, 0, 0), (R, 0, T), (0, 0, T). Notice that Eq. (64) 
is gauge invariant and to further proceed we must fix the gauge. Under the Feynmann gauge [5S], it becomes 



W h [C] cx J D[A h \ex$ \ 



9 ( .2 (' / /'-^') '",.'V 



(Bl) 



Since we are interested physics at large scales, we pass to the continuum limit, obtaining the gauge field propagator 
read 



dud 2 k 



ii i\ f dud 

rWx-x,T-r) = y ^ 



p 2 & 

,ik-(x— x )— iu(r— t') 



1 



|k| 2 4tt ^|x-x'| 2 + \t-t'\ 2 ' 



(B2) 



Then, because of T 3> R diagrams with a propagator line starting from the rectangular sides in the spatial direction 
are penalized by (some power of) a small factor JJ/T « 1 and thereby are negligible. As a result, to the leading 
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order expansion in e\ Eq. (Bl I is dominated by the two diagrams shown in Fig. [4] By ignoring the irrelevant overall 
numerical factor, which is the same for both diagrams, they give 

diagram (a) = -2 x I J D[A h ] £ drdr'A^O, 0, t)A%{0, 0, r') exp £ (d^)* 

= -/ dTdr'D 0O (0,0,T-T') = -^- [ drdT 1 —^— 
Jo 47r Jo \t-t'\ 



re; 
' 2tt 



2 r i 



dx In 



xT 



(B3) 







with t* an ultraviolet cutoff in the imaginary time direction, and 



diagram (b) = 2 x 1 1 jf drdr'A h o (0, 0, r)i4§(ij, 0, r') exp £ (d M Aj) 

= [ T dTdT'D oo (R,0,T~T') = ^ [ T dTdr'—=^= 
Jo 47r + I 7 " - T I 



o2 /" T 



4-7T 



dr In 



T-r+ ^(r-^^i? 2 Te 2 f 1 x+^/x 2 + {R/ry 



-t + Vt 2 + R 2 



in 



dx In ■ 



+ v/x 2 + [R/Tf 



t^oo Te\ f 1 , , 2Tx 

~ / da; m , 

2tt J R 



(B4) 



respectively. Notice that in the first line of Eqs. (B3) and (B4), the factor 2 is the combinatorial factor. Adding these 
two terms together, we find 



diagram (a) + diagram (b) 



t^oo Te? , R Te 2 , R 

~ m = In — . 

2tt 2t* 2vr ft* 



(B5) 



which is Eq. (65 1, with ft^ = 2r* an ultraviolet cutoff in the spatial direction. Notice that these two diagrams suffer 
logarithmical divergences in T , which, however, cancel each other exactly upon adding them together. 



(a) 



(b) 





> — D uv (x x ) 



x x ^ v 

FIG. 4: Diagrams contributing to the Wilson loop for T — > 00. 
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Appendix C: Spinon superfluids 



From Eqs. (42) and (70) we obtain 



(CI) 



If we re-arrange the components of the vector as follows: = (^oi, ^02, *03, *04) T — > (*oi) ^04) ^02, *03) T > 
and further introduce two two-component vectors, z a , a =T, i, defined in the sector introduced by the sublattice 
structure, then ^0 — \ Z J) T ■ (This two-component structure is defined in the spin sector.) Consequently, we may 
rewrite Eq. (ICl]) as 



X s 

-2J S 

X s 
-2J, 



Solving these two equations we obtain 



-2J S 
X s 

-2J S 
X s 



Zf + u 2 {z\z^ + z^z^z^ = , 
H + u 2 {z\z^ + z\z i )z i = . 

n t +n i = (2J S - A s )/(2u 2 ) , 



(C2) 
(C3) 

(C4) 



where n a and a are homogeneous in spacetime. This saddle point solution gives z\Tj^z a = with £3 = diag(l, — 1) 
defined on the sector introduced by the sublattice structure. As a result, we find that the total spin polarization 
vanishes namely Eq. ( 71 1 . 
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